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Abstract 

Inflation of a universe filled by an elastic continuous medium is considered. Elastic 
inflation, while capable of describing current observations, is qualitatively different from 
standard inflationary models. The scalar and tensor modes keep evolving after crossing the 
horizon. Due to this superhorizon evolution, the amplitude of the tensor mode (inflationary 
gravitational wave) is not simply proportional to the energy density of the universe at the 
time of horizon crossing. The spectral index of tensor modes can be positive. 
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I. INTRODUCTION 



The cosmo logical theory, and cosmological perturbation theory in particular [1, 2], have 
been confirmed observationally only when applied to an ideal fluid. Yet inflation and infla- 
tionary perturbations are calculated assuming that the universe was dominated by a scalar 
field [2, 3, 4], or a higher-derivative gravity [5], a scalar field with non-minimal kinetic term 
[6], a higher-derivative scalar field [7], a self-coupled vector field [8], or by a self-coupled 
three form [9]. 

Since physics at inflationary energies is unknown, it makes sense to consider a fluid-like 
description of inflation, which may ultimately turn out to be closer to reality than any of the 
above-listed field models. An ideal fluid cannot drive inflation, but an ideal elastic medium 
can. Elastic medium was already considered as a model of dark energy [10]. 

We describe the relativistic elastic medium in §2, and the elastic inflation in §3, where 
we also introduce the speeds of longitudinal and transverse sound waves. We calculate the 
inflationary tensor spectrum in §4, scalar spectrum in §5, and discuss the results in §6. 

II. ELASTIC MEDIUM 

Continuous elastic medium is a collection of time-like world lines filling the space-time 
("congruence", see [10] and references therein). The action of elastic medium is given by 
a time integral of potential energy along each of the world lines of the congruence. The 
potential is an arbitrary function of spatial distances between neighboring world lines. 

It is convenient, at least for our purposes, to choose spatial coordinates coincident with 
the markers of the world lines, so that the world lines of the congruence become the lines 
of constant x\ % = 1,2,3. Then the action of the elastic medium can be written in the 
following form 




(1) 



Here V is an arbitrary function of the invariants of the spatial metric 



lij — 9ij + 
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(2) 
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As the three invariants, one may take 

h = la, t 2 = ifijlju 7 = det(7ij). (3) 

We will not go beyond the first-order perturbation theory (will not calculate the non- 
Gaussianities for elastic inflation). Then, given that deviations from isotropy are small 
and 7ij is approximately diagonal, to sufficient accuracy, we may express the third invariant 
as a function of the other two. We choose the traces t\ and t 2 as the variables of elastic 
energy, V = V(t 1 ,t 2 ). All the formulas written below are valid only up to first order in 
metric perturbation. 

The energy-momentum tensor of elastic medium is obtained by varying the action with 
respect to g^ v . (The non-covariant form of the following expressions is due to the special 
choice of coordinates. All our final results will be expressed in a coordinate-free language.) 

T oo = Vg oo^ (4) 
T 0i = Vg 0i + 2V 1 ^-AV 2 ^^, (5) 

T ik = y g ik _ 2Vi5 . k + A y 2g . k (6) 

Here and in the following, V\ = d tl V, V 2 = d t2 V. With 87rG = 1, the Einstein equations are 

Gt = T$. (7) 

When the energy-momentum tensor is known, one calculates the unperturbed inflation and 
inflationary perturbations along standard lines [2]. 

III. ELASTIC INFLATION 

The background metric is flat FRW 

ds 2 = a 2 (r]) (dif — 5ijdx % dx^ , 
and the corresponding Friedmann equations are 

3a- 2 H 2 = V =e, 
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(8) 
(9) 



a~ 2 (2H' + H 2 ) = V + 2a 2 V l + 4aV 2 = -p, (10) 

here and in what follows H = a'/a and the prime denotes d v . We have also defined the 
unperturbed energy density and pressure. 

Inflationary stage requires p ~ — e. On the other hand, we will show that 

^ = c 2 --c 2 (11) 
dt s 3 1 J 

where c s is the speed of scalar (longitudinal) sound and c v is the speed of vector (transverse) 
sound. Thus, inflation requires 

c 2 -^c>-l, 0<c 2 <l, 0<c 2 <l. (12) 

These requirements constrain inflationary potentials V(ti,t 2 ). 

In the above equations, c s and c v are the short- wavelength limits of the propagation speeds 
of scalar and vector perturbations. These can be calculated neglecting gravity (assuming a 
pure gauge metric and neglecting the expansion of the universe). As shown in Appendix A, 

2 _ 3Fi + Wa 2 V 2 + 2a 2 V n + 8a 4 V 12 + 8a 6 V 22 

° s ~ V 1 + 2aW 2 ' {l6) 

c 2 = - (14) 
V 1 + 2a 2 V 2 ' { } 

where Vn,...are second derivatives. Equation (11) follows. 

For c v = 0, one can show that V = V{^). Then the energy depends only on the change of 
volume, and the elastic model reduces to an ideal fluid. We will assume that elastic inflation 
ends when c v vanishes. If at this moment also c s = 0, then one needs to add reheating. But 
reheating is not necessary for elastic inflation, one may assume that at the end of inflation 
c 2 s — 1/3 - inflation ends when an ultrarelativistic fluid loses elasticity. 

A simple elastic inflation model, which we use for numerical examples, assumes constant 
sound speeds c S j, c V i during the inflationary epoch, satisfying c 2 si — (4/3)c^ = —1. Inflation 
ends at the energy density e r . After inflation, for e < e r , we also assume constant sound 
speeds, c sr , and c vr = 0. Equation (11) gives the pressure 

p = c 2 sr e, e < e r , p = (1 + c 2 sr )e r — e, e > e r . (15) 
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Friedmann equations give energy 



e = e r + 3(1 + c 2 sr )e r N, (16) 

as a function of the number of e-foldings until the end of inflation, N = — In a + const. 
We will see that this artificial "constant speeds" model satisfies all current observational 
constraints on inflation. 

IV. TENSOR MODES 

We start with tensor modes, because the calculation of the spectrum is much simpler 
than for scalars. The perturbed metric is 

ds 2 = a 2 (drf - (5ij + h ij )dx i dxP > ) , (17) 

where is a symmetric pure tensor, with ha = and djhy = 0. Equations (4-6) give the 
energy-momentum tensor perturbation 

5T] = 2a 2 (V 1 + 4a 2 V 2 )h lj . (18) 

The perturbation of the Einstein tensor is 

5G) = - X -ar 2 {b!' l3 + mh! l3 + k 2 hij) , (19) 

where k is the comoving wavenumber. Dropping the indices, we get the following tensor 
mode equation 

h" + 2Hh' + k 2 h + 2c 2 v a 2 (e + p)h = 0. (20) 

It is convenient to write the mode equation using cosmic time, dt = adrj. Denoting the time 
derivative by the dot and introducing the Hubble constant H = a/a, we have 

k 2 

h + 3Hh+— h + 2c 2 v (e + p)h = 0. (21) 
a z 

For zero elasticity, c v = 0, the above equation reduces to the standard one [2]. 

Denote by Fx(k) the elastic transfer function for tensor modes defined as the large time 
limit of the ratio of the amplitudes of the mode described by equation (21) and the mode 
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described by (21) with the same a(t) but with c v = 0. The late-time ratio of amplitudes is 
calculated for the modes with equal early-time (Ha <C k) amplitudes. 

The transfer function can be calculated analytically assuming slow-roll inflation. The 
effect of elasticity remains small before horizon crossing and leads to a slow roll of h after 
horizon crossing, according to 

3Hh + 2c 2 v (e + p)h = 0. (22) 

The last equation gives h oc F T (k) with 

Fr(k) ~ exp(-^£ fe ^) (23) 

Here e r is the energy density at the end of inflation, e& is the energy density at the time 
when the k mode crosses horizon, that is when Ha = k. 

The slow-roll condition c 2 (e + p) <£ will be violated near the end of inflation, leading to 
a multiplicative k-independent error of order few in expression (23). The correct expression 
for the transfer function is 

F T (k) = F T exp (- \ £ ^ ) (24) 

where F T is a constant factor of order few. 

The tensor perturbation spectrum is (here and for the scalar spectrum we use the defini- 
tions used in [11]) 

>h,u\ 128 e k / 4 p c 2 v de N 



p, w = -^s; exp "si i-J (25) 



The tensor spectral index (the tilt) is 



For zero elasticity, c v = 0, this reduces to the standard expression [11]. However, elastic 
inflation requires c 2 s — —1, giving 

n T = 3c 2 s (l + . (27) 

For elastic inflation the tensor index can be positive. 
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For the constant speeds model of elastic inflation described at the end of §3, we get the 
tensor power 

and the tilt 

n T = §- (29) 

V. SCALAR MODES 

Generic scalar perturbations of the metric are, [2], 

ds 2 = a 2 ( (1 + 2 < p)drj 2 - 2d i Bdr ] dx i - - 2^% + 2d i d j E)dx i dx j ) . (30) 

For this metric one finds the perturbed energy-momentum tensor from (4-6) and the per- 
turbed Einstein tensor from [2]. We do this in Appendix B, here we just give a summary of 
results. Define a new variable ( by 

- (l + 4^) C = * + ^(*' + (l- |<i)H.), (31) 

where ^ = ip — Ti(B — E') is the gauge-invariant potential. After inflation, c v = 0, ( 
becomes the standard curvature perturbation (on superhorizon scales, for c v = 0, ( is the 
perturbation of the scale factor on uniform energy hypersurfaces). Our ( also coincides 
with the usual ( at very early time when k ^> H, allowing us to use standard results when 
calculating quantum fluctuations. 
The ^-perturbation satisfies 

C + -C + -VC + < = o. (32) 



z a 2 



Here 



_ a 3 (e + p) 

z = -&p- (33) 



coincides with the corresponding function for an ideal fluid case [11]. The other function 



w = Ac 2 v H(2^-2^-l 6 -) (34) 



c v c s 2 e 
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turns to zero in the ideal fluid limit. 

The scalar perturbation spectrum is calculated from (32) similarly to the tensor mode 
case. Using [11], we find 

P<(k) = '4 Fl , 1 } , exp{- t r ^ ( 1 + 4^ - 4 d -p^\ }. (35 ) 

v ; 9 5 M P c s {l+p k /e k ) P1 3Je r e \ dine dine) 1 v ; 

The scalar tilt is 



— U + + (36) 



^(*) = ^J^ — (-)'" (37) 



For the constant speeds inflation described at the end of §3, the scalar power is 

16 2 3N, 
9 s M P c sl 
and the tilt is 

ns= - (38) 

VI. DISCUSSION 

Elastic Inflation is capable of describing available observations. Currently, there exist 
three inflationary observables [12, 13]: the tensor to scalar ratio, the scalar tilt, and non- 
Gaussianity (all upper bounds). The scalar power, which has been measured to a great 
accuracy, simply gives the energy scale of inflation, and does not constrain phenomenological 
models like ours. 

We have not calculated the non-Gaussianity, but one expects it to be either ~ 1/N ([14]) 
or ~ 1 ([15]), easily passing the observational constraints. The other two constraints are 
satisfied even by the artificial constant speeds model of elastic inflation described at the end 
of §3. This model gives a negative scalar tilt —(1 — c 2 si )/N ', which for N ~ 60 is consistent 
with observations for any < c si < 1. The tensor to scalar ratio ~ Sc si /N is also consistent 
with current observations for arbitrary < c si < 1. 

Summarizing: 

(i) Elastic Inflation is described by equations (9, 10, 11, 12). 

(ii) The tensor and scalar spectra generated by elastic inflation are given by equations 
(25, 26, 35, 36), where the constant factors of order unity Ft,s should be calculated by 
numerical integration of the mode equations (21, 32). 
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APPENDIX A: SOUND SPEEDS 

We consider perturbations on Minkowski background ds 2 = dt 2 — 5ijdx l dx j . For short- 
wavelength perturbations we neglect gravity, meaning that the metric perturbation is a pure 
gauge, x % — > x % — £ l , giving the perturbed metric 

fl'oo = 1, 9oi = 9ij = Sij - di£j - dj£i. (Al) 

The perturbed traces and 4-Jacobian are 

ti = 3 + 20 fc &, *2 = 3 + 40 fc &, V=g = l + d k £ k . (A2) 

The perturbed energy-momentum tensor is calculated from (4-6): 

T° = V + 5V, Tf = (2V 1 + AV 2 )i, T* = -2(V + 2V 1 + AV 2 )i, (A3) 

T k = (5V + 25V, + A5V 2 )5 lk + (2V, + 8V 2 )(d^ k - <9 fc &), (A4) 

where V\ = d tl V,... and the perturbations of potential energy and its derivatives are given 
by 

SV = (2V, + W 2 )d k &, 5V = (2V„ + W 12 )d k &, 5V 2 = (2V 21 + W 22 )d k £ k . (A5) 

The equation of motion 

= = -L^d^y/^TS) - \(d v g a p)T^ (A6) 

V 9 z 

gives 

(2V, + AV 2 % + d k (5V + 25V, + A5V 2 ) + (2V, + W 2 )(d% + d k d n £ n ). (A7) 

For vector modes (d n ^ n = 0) we get the propagation speed (14), with a=l. For the scalar 
mode (£ n = d n x), we get the speed (13). 
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APPENDIX B: SCALAR PERTURBATIONS 

The perturbed metric is 

ds 2 = a 2 ( (1 + 2<p)dri 2 - 2d i Bdr ] dx i - - + 2d i d j E)dx i dx j ) . (Bl) 

Perturbation of the energy-momentum tensor is found from (4-6) 

5T ° = 5V, 57? = (2aVi + 4aV 2 )^5, (B2) 

5T k = (5V + 2a 2 5V 1 + Aa i 5V 2 )5 lk + (4a V x + 16aV 2 )(^9 fc E - ^5 ik ). (B3) 
Here the perturbations of potential energy and its derivatives are 

SV = V 1 8t 1 + V 2 5t 2 , 5V 1 = V n St 1 + V 12 5t 2 , 5V 2 = V 21 5h + V 22 5t 2 , (B4) 

and the perturbations of the traces are 

5ti = -2a 2 (3ip + k 2 E), 5t 2 = -4a 4 (3^ + k 2 E), (B5) 

where k is the comoving wavenumber. 

With Einstein tensor written in [2], the Einstein equations are 

- m 2 <p - SHip' - k 2 ij + k 2 H{B - E') = \a 2 5V, (B6) 

H<j) + ^' = a 4 (V 1 + 2a 2 V 2 )B, (B7) 

k 2 

(2H' + H 2 )(f) + Htf + V" + - —D = 

Z 

-\a 2 {bV + 2a 2 5V l + Aa 4 5V 2 ) + 2a 4 (Vi + 4aV 2 )^, (B8) 
z 

D = 4a 4 (Vi + 4a 2 V 2 )E, (B9) 

D = <f> - i/j + (B - E')' + 2H(B-E'). (BIO) 

This system is redundant. Due to time reparametrization invariance one equation can be 
dropped. Equivalently, one gauge condition may be imposed. We set 

B = 0. (Bll) 

10 



Then (B7) gives 

ft0 + -0' = O. (B12) 
Introducing a gauge invariant potential (see [2] for discussion) 

^ = vb + HE\ (B13) 

we get a system of two first-order ODEs 

+ H^ = a 4 (Vi + 2a 2 V 2 )E' - Aa^^ + Aa 2 V 2 )HE (B14) 

k 2 ^ = cfiVx + 2aV 2 )(3* - 3HE' + k 2 E). (B15) 

Now we want to get an equation for the curvature perturbation (. More precisely, we 
would like to define a variable which becomes curvature perturbation in the ideal fluid limit 
c v = (because this is what we want to calculate), and also coincides with the standard ( 
variable at early times (this will allow to use standard results to normalize the mode). This 
can be done in the following (rather tedious) way. 

One notes that the system (B14, B15) can be written as 

(2a* + a 3 (e + p)E)' = -3c 2 Ha 3 (e + p)E, (B16) 

k 2 (2a^ + a 3 (e + p)£) = -3a 3 (e + - HE'). (B17) 
With the help of [11], one defines a possible variable (: 

( = ^( E + ^r--^)- (bis) 

6 a z {e + p) 
The definitions and equations given in the main text follow. 
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